subroutine loop1(umaru, umc, ump, aru, chic, chip, mza, mua) 

implicit none

real(kind=8), intent(inout) :: umaru  (mza,mua)
real(kind=4), intent(in) :: umc  (mza,mua)
real(kind=4), intent(in) :: ump  (mza,mua)
real(kind=4), intent(in) :: aru  (mza,mua)
real(kind=4), intent(in) :: chic, chip 
integer, intent(in) :: mza, mua
integer :: iu,k

umaru=0.0

do iu = 2,mua
   do k = 1,mza-1
      !umaru(k,iu) = (chic * umc(k,iu) + chip * ump(k,iu)) * aru(k,iu)
      umaru(k,iu) = (umc(k,iu) + ump(k,iu)) * aru(k,iu)
   enddo
enddo

end subroutine loop1
